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Abstract 



The paper studies an implementation methodology for partial and 
disjunctive stable models where partiality and disjunctions are un- 
folded from a logic program so that an implementation of stable models 
for normal (disjunction-free) programs can be used as the core infer- 
ence engine. The unfolding is done in two separate steps. Firstly, it is 
shown that partial stable models can be captured by total stable mod- 
els using a simple linear and modular program transformation. Hence, 
reasoning tasks concerning partial stable models can be solved using 
an implementation of total stable models. Disjunctive partial stable 
models have been lacking implementations which now become avail- 
able as the translation handles also the disjunctive case. Secondly, it 
is shown how total stable models of disjunctive programs can be de- 
termined by computing stable models for normal programs. Hence, an 
implementation of stable models of normal programs can be used as a 
core engine for implementing disjunctive programs. The feasibility of 
the approach is demonstrated by constructing a system for comput- 
ing stable models of disjunctive programs using the SMODELS system 
as the core engine. The performance of the resulting system is com- 
pared to that of DLV which is a state-of-the-art system for disjunctive 
programs. 

1 INTRODUCTION 

Implementation techniques for declarative semantics of logic programs have 
advanced considerably during the last years. For example, the XSB sys- 
tem [40j is a WAM-based full logic programming system supporting the well- 
founded semantics. In addition to this kind of a skeptical approach that is 
based on query evaluation also a credulous approach focusing on computing 
models of logic programs is gaining popularity. This work has been centered 
around the stable model semantics [13 HE]- There are reasonably efficient 
implementations available for computing stable models for disjunctive and 
normal (disjunction- free) programs, e.g., DLV |24], SMODELS [ifJIlEl, CMOD- 
ELS [T], and ASS AT [2^1. The implementations have provided a basis for a 
new paradigm for logic programming called answer set programming (a term 
coined by Vladimir Lifschitz). The basic idea is that a problem is solved by 
devising a logic program such that the stable models of the program provide 
the answers to the problem, i.e., solving the problem is reduced to a stable 
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model computation task [23 E2 EH 1131 El ■ This approach has led to interest- 
ing applications in areas such as planning [HI HUE], model checking [3 01 I2U]. 
and software configuration [19] . 

This paper addresses two issues in the stable model semantics: partial- 
ity and disjunctions. The idea is to develop methodology such that efficient 
procedures for computing (total) stable models that are emerging can be ex- 
ploited when dealing with partial stable models and disjunctive programs. 
Sometimes it is natural to use partial stable models to represent a domain. 
Even when working with total stable models, partial stable models could 
be useful, e.g., for debugging purposes to show what is wrong in a program 
without any total stable models. However, little has been done on imple- 
menting the computation of partial stable models and most of the work has 
focused on query evaluation w.r.t. the well-founded semantics. In the paper 
we show that total stable models can capture partial stable models using 
a simple linear program transformation. This transformation works also in 
the disjunctive case showing that implementations of total stable models, 
e.g. dlv, can be used for computing partial stable models. Using a suitable 
transformation of queries, a mechanism for query answering can be realized 
as well. 

Our translation is interesting in many respects. First, it should be noted 
that the translation does not follow directly from the complexity results al- 
ready available. It has been shown, e.g., that the problem of deciding whether 
a query is contained in some model (possibility inference) is Sf-complete for 
both partial and total stable models [l2l[TH]. This implies that there exists a 
polynomial time reduction from possibility inference w.r.t. partial models to 
possibility inference w.r.t. total models. However, this kind of a translation 
is guaranteed to preserve only the yes/no answer to the possibility inference 
problem. Second, not all translations are satisfactory from a computational 
point of view. In practice, when a program is compiled into another form 
to be executed, certain computational properties of the translation play an 
important role: 

• efficiency of the compilation (in which order of polynomial), 

• modularity (are independent, separate compilations of parts of a pro- 
gram possible), and 

• structural preservation (are the composition and intuition of the origi- 
nal program preserved so that debugging and understanding of runtime 
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behavior are made possible). 

All this points to the importance of finding good translation methods to 
enable the use of an existing inference engine to solve other interesting prob- 
lems. 

The efficiency of procedures for computing stable models of normal pro- 
grams has increased substantially in recent years. An interesting possibility 
to exploit the computational power of such a procedure is to use it as a core 
engine for implementing other reasoning systems. In this paper, we follow 
this approach and develop a method for reducing stable model computation 
of disjunctive programs to the problem of determining stable models for nor- 
mal programs. This is non-trivial as deciding whether a disjunctive program 
has a stable model is S^-complete [l2j whereas the problem is NP-complete 
in the non-disjunctive case [Hi]. The method has been implemented using 
the SMODELS system |46l 145) as the core engine. The performance of the im- 
plementation is compared to that of DLV, which is a state-of-the-art system 
for computing stable models for disjunctive programs. 

There are a number of novelties in the work. Maximal partial stable 
models for normal programs are known as regular models, M-stable models, 
and preferred extensions [El EH1 150) . Although this semantics has a sound 
and complete top-down query answering procedure [I0J HOJ EH], so far very 
little effort has been given to a serious implementation. For disjunctive pro- 
grams, to our knowledge, no implementation has ever been attempted. As 
a result, we obtain (perhaps) the first scalable implementation of the regu- 
lar model/preferred extension semantics, and the first implementation ever 
for partial stable model semantics for disjunctive programs. Our techni- 
cal work on the relationship between stable and partial stable models via a 
translational approach provides a compelling argument for the naturalness 
of partial stable models: stable models and partial stable models share the 
same notion of unfoundedness, carefully studied earlier in [I4j[26j. Finally, 
we demonstrate how key tasks in computing disjunctive stable models can 
be reduced to stable model computation for normal programs by suitable 
program transformations. In particular, we develop techniques for mapping 
a disjunctive program into a normal one such that the set of stable models 
of the normal program covers the set of stable models of the disjunctive one 
and in many case even coincides with it. Moreover, we devise a method 
where the stability of a model candidate for a disjunctive program can be 
determined by transforming the disjunctive program into a normal one and 
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checking the existence of a stable model for it. Finally, in the experimental 
part of this paper, we present a new way of encoding quantified Boolean 
formulas as disjunctive logic programs. This transformation is more econom- 
ical in the number of propositional atoms and disjunctive rules than earlier 
transformations presented in the literature [T2*l 125]. 

The rest of the paper is structured as follows. We first review the basic 
definitions and concepts in Section El It is then shown in Section El that par- 
tial stable models can be captured with total stable models using a simple 
program transformation. In Section El we describe the method for com- 
puting disjunctive stable models using an implementation of non-disjunctive 
programs as a core engine. After this, we present some experimental results 
in Section El and finish with concluding remarks in Section H3 

As a comment on the historical development of the translation given in 
Section 3, the characterization of partial stable models as stable models of 
the transformed program was first sketched for normal programs in a proof 
by Schlipf (42l Theorem 3.2]. For disjunctive programs, it was discovered and 
proven in [IB], and independently in [21]. In the current paper we present a 
proof based on unfounded sets, which was given in [21], as this proof reveals 
some of the properties of unfounded sets which are of interest in their own 
right. Yet another approach to computing the partial stable models of a 
disjunctive program based on a program transformation has been developed 
by Ruiz and Minker [HI|: a disjunctive program P is translated into a positive 
disjunctive program P 3S with constraints, the 3 S -transformation of P, such 
that the total minimal models of P 3S that additionally fulfill the constraints 
coincide with the partial stable models of P. 

2 DEFINITIONS AND NOTATIONS 

A disjunctive logic program P (or, just disjunctive program P) is a set of rules 
of the form 

ai V • • • V a k <— h, . . . , b m , ~ci, . . . , ~c n (1) 

where k > 1, m, n > and a^'s, fej's and c,'s are atoms from the Herbrand 
base Hb(P) 1 of P. Let us also distinguish subclasses of disjunctive programs. 

1 For the sake of convenience, we assume that a given program P is already instantiated 
by the underlying Herbrand universe, and is thus ground. 
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If k = 1 for each rule of P, then P is a disjunction-free or normal program. 
If n = for each rule of P, then P is called positive. 

Literals are either atoms from Hb(P) or expressions of the form ~a where 
a G Hb(P). For a set of atoms A C Hb(P), we define ~A as {~a | a G A}. 
Let us introduce a shorthand A <— B, ~C for rules where A ^ 0, B and C 
are subsets of Hb(P). In harmony with (JIJ, the set of atoms A in the head 
of the rule is interpreted disjunctively while the set of literals B U ~C in 
the body of the rule is interpreted conjunctively. We wish to further simplify 
the notation A <— B, ~C in some particular cases. When A, B or C is a 
singleton {a}, we write a instead of {a}. If B = or C = we omit P and 
~C (respectively) as well as the separating comma in the body of the rule. 

2.1 PARTIAL AND TOTAL MODELS 

We review the basic model-theoretic concepts by following the presentation 
in [11]. Let P be any disjunctive program. A partial interpretation I for P is 
a pair (T, F) of subsets of Hb(P) such that T R P = 0. The atoms in the sets 
J* = T, J f = P and J u = Hb(P) — (TUP) are considered to be true, false, 
and undefined, respectively. We introduce constants t, f, and u, to denote 
the respective three truth values. A partial interpretation I for P is a total 
interpretation for P whenever I u = 0, i.e., if every atom of Hb(P) is either 
true or false. When no confusion arises, we use P alone to specify a total 
interpretation / for P (then P" = Hb(P) — P and I u = hold). 

Given a partial interpretation for P, the truth values of atoms are de- 
termined by P, I f and I u as explained above while t, f and u have their 
fixed truth values. For more complex logical expressions P, we use 1(E) to 
denote the truth value of E in J. The value J(~a) is defined to be t, f, 
or u whenever 1(a) is f, t, or u, respectively. To handle conjunctions and 
disjunctions, we introduce an ordering on the three truth values by setting 
f < u < t. By default, a set of literals L = {h, . . . , l n } denotes the conjunc- 
tion li A • • • A l n while \J L denotes the corresponding disjunction l\ V • • • V l n . 
The truth values I(L) and /(V L) are defined as the respective minimum 
and maximum among the truth values 1(h), ■ ■ ■ , 1(h)- A rule A <— B, ~C 
is satisfied in I if and only if I(\f A) > I(B U ~C). A partial interpretation 
M for P is a partial model of P if all rules of P are satisfied in M, and for 
a total model, also M u = holds. Let us then introduce an ordering among 
partial models of a disjunctive program: Mi < M 2 if and only if Mi* C M 2 * 
and Mi f D M 2 f . A partial model M of P is a minimal one if there is no 
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partial model M' of P such that M' < M (i.e., M' < M and M' ^ M). In 
case of total models, we have N\ < N 2 if and only if Ni C N 2 . Moreover, 
a total model N of P is considered to be a minimal one if there is no total 
model N' of P such that N' C N. 

2.2 STABLE MODELS 

Given a partial interpretation / for a disjunctive program P, we define a 
reduction of P as follows: 

P 1 = {A <- B | A <- 5, ~C G P and C C J f }. 

Note that this transformation coincides with the Gelfond-Lifschitz reduction 
of P (the GL-reduction of P) when / is a total interpretation. 

Definition 2.1 (Total stable model) A total interpretation N for a dis- 
junctive program P is a stable model if and only if N is a minimal total model 
ofP N . ' 

The original definition of partial stable models |3E1 E5H is based on a 
weaker reduction. Given a disjunctive program P and an interpretation 
/, the reduction Pi is the set of rules obtained from P by replacing any 
~c in the body of a rule by /(~c). As noted in [SHI, the only practical 
difference between P 1 and Pi is that Pj has rules that correspond to rules 
of A <- B, ~C G P satisfying J(~C) = u. Note that if J(~C) = t, then 
A <— B G P 1 , and if J(~C) = f, then the partial models of Pi are not 
constrained by the rule included in Pj. 

Definition 2.2 (Partial stable model) A partial interpretation M for a 
disjunctive program P is a partial stable model of P if and only if M is a 
minimal partial model of Pm- 

In the above definition, the relation between M and Pm is similar to 
the one for the total stable model, both for the purpose of preserving the 
stability condition. While maximizing falsity and minimizing true atoms, a 
partial stable model does not insist that every atom must be either true or 
false. (Partial) stable models are intimately related to unfounded sets [THI2B] . 

Definition 2.3 (Unfounded sets) Let I be a partial interpretation for a 
disjunctive program P . A set U C Hb(P) of ground atoms is an unfounded 
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set for P w.r.t. I , if at least one of the following conditions holds for each 
rule A <- B,~C G P such that AnU 7^ 0: 

UF1: Bnl f ^ orCH/V 0, 

UF2: Bn U ^ 0, or 

UF3: (A-U)H (P U J u ) ^ 0. 

^4n unfounded set U for P w.r.t. I is I -consistent if and only if U fl P — 0. 

The conditions UF1 and UF3 above coincide with the conditions 

I{B U ~C) = f and J(V (-4 - 17)) ^ f , 

respectively. The intuition is that the atoms of an unfounded set U can be 
assumed to be false without violating the satisfiability of any rule A <— B>, ~C 
of the program whose head contains some atoms of U. For any such rule, 
either the rule body is false in I (UF1), or the rule body can be falsified 
by falsifying the atoms in U (UF2), or the head of the rule is not false in I 
(UF3). In particular, unfounded sets w.r.t. partial/total models can be used 
for constructing smaller partial/total models (recall the definition of minimal 
partial and total models) in a way that is made precise by what follows. 

Lemma 2.4 Let M = (T, F) be a partial model of a positive disjunctive 
program P and U an unfounded set for P w.r.t. M . Then, if M is total or 
U is M-consistent, M' = (T — U , F U U) is a partial model of P. 

PROOF. Let M, P, U and M' be defined as above. Additionally, we assume 
that (a) M is total, or (b) U is M-consistent. Let us then assume that some 
rule A <- B of P is not satisfied in M' which means that M'(\J A) < M'(B). 
Thus (i) M'(V A) < t and M'(B) = t, or (ii) M'{\j A) = f and M'(B) = u. 
Our proof splits in two separate threads. 

I. Assume that AnU = holds. Consider the case (i). Now M'(\f A) < t 
implies A n (T - U) =0. Since A n U = holds, too, we obtain 
A n T = so that M(\j A) < t. On the other hand, M'(B) = t 
implies B C T — U. Thus B C T and M(B) = t holds as well. 
But then M(\J A) < M(B), a contradiction. The case (ii) is analyzed 
next. Now M'(\J A) = f implies A C F U U as well as A C F, since 
AnU = 0. Thus M(V A) = f. Moreover, from M'(B) = u we obtain 
Bn(FUU) = 0. Thus we obtain BnF ^ so that M{B) > f holds. To 
conclude, we have established that M(\J A) < M(B), a contradiction. 
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II. Otherwise A fl U ^ holds. Then at least one of the unfoundedness 
conditions is applicable to A <— B. If UF1 is, M(B) = f holds. It 
follows that B C F and B C F U [/. Thus M'(5) = f contradicting 
both (i) and (ii). If UF2 is applicable, we have B nU ^ 0. It follows 
that B n (F U 17) 7^ so that Af'(S) = f , a contradiction. 

Thus UF3 must apply, i.e., M{\J{A-U)) > f holds. Let us then 
consider cases (a) and (b) separately. 

(a) If M is total, we have necessarily M(\f (A — U)) = t. This implies 
that some atom a e A — U belongs to T. Thus also a G T — U 
and M'(\J A) = t, a contradiction with both (i) and (ii). 

(b) If U is M-consistent, we have U n T = 0. By M(V (A -£/))= t 
there is an atom a & A — U such that a ^ F. Then a ^ F U U 
which implies M f (\f A) > f. Thus (ii) is impossible and (i) implies 
M'{\JA) = u and M'{B) = t. It follows that A(l(T — U) = 
and B C T — U. Since C/ fl T = 0, the former implies A fl T = 
while the latter implies that B CT. Consequently, M(\f A) < t 
and M(77) = t, i.e., M(\/ A) < M(B), a contradiction. 

□ 

Let us yet emphasize the content of Lemma [2~H when M is total (and [/ 
need not be M-consistent). Then M' = M — U is also a total model of P. A 
couple of examples on unfounded sets follow. 

Example 2.5 Consider a disjunctive program 

P = { a V b <- c, ~a} 

and an interpretation I = (0, {a}). The only rule in P has its body undefined 
in J, hence UF1 is not applicable. The set {a} is unfounded w.r.t. I since b 
is undefined in I and not in the set, hence UF3 is applicable. On the other 
hand, the set {b} is not unfounded w.r.t. I whereas {c} is unfounded w.r.t. 
7. Once c belongs to an unfounded set, the atoms a and b can both get in 
due to UF2. Hence, we have U = {a, b, c} as an unfounded set w.r.t. 7. 

Comparing U with 7 f , we find that 7 f does not maximize the atoms that 
should be false. This program P has exactly one stable model (which is also 
a partial stable model) in which all three atoms are false. □ 
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Unlike the case for normal programs, the union of unfounded sets may 
not be an unfounded set. 

Example 2.6 Consider a program P containing only one rule 

a V b <- 

and an interpretation / = ({a, 6},0). The program has two non-empty un- 
founded sets w.r.t. J, {a} and {b}. Either a or b depends on the other one not 
in the set for UF3 to be applicable. However, UF3 becomes not applicable 
when both a and b are in, thus the union {a, b} is not an unfounded set. □ 

An interpretation / for P becomes particularly interesting when the union 
of all unfounded sets U for P w.r.t. I is also an unfounded set for P w.r.t. 
I. In this case, the program P possesses the greatest unfounded set U for P 
w.r.t. I. 

Definition 2.7 A total interpretation I is said to be unfounded free for a 
program P if and only if there is no unfounded set U for P w.r.t. I such that 

(/n/VI. 

The notion of unfounded freeness captures the stable model beautifully. 

Theorem 2.8 |26j Let M be a total interpretation for a disjunctive program 
P. Then, the following are equivalent 

• M is a stable model of P. 

• M f is the greatest unfounded set for P w.r.t. M. 

• M is unfounded free for P. 

On the other hand, Eiter et al. [TJ] show that partial stable models can 
be defined essentially without reference to three-valued logic. 

Theorem 2.9 [14| If M is a partial interpretation for a disjunctive program 
P, then M is a partial stable model of P if and only if 

• M t is a minimal total model of P M and 

• M f is a maximal M -consistent unfounded set for P w.r.t. M. 
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The first condition in the theorem is called foundedness in [H]. The dif- 
ferences between these two theorems are quite subtle. The strictness of stable 
models enforces a simpler relationship between stable models and unfounded 
sets. Therefore, neither maximality nor consistency nor foundedness need 
be explicitly stated. The characterization of partial stable models in Theo- 
rem EH accounts for a more reflexible situation: since M may not be a total 
model, maximality should extend the set of false atoms as much as possible 
without causing inconsistency. However, maximality and consistency are still 
not strong enough. 

Example 2.10 Consider a disjunctive program 

P = { a V b «- ~a} 

and an interpretation / = ({a}, {&}) which is total so that the definition of 
unfoundedness makes no difference in valuation under /. Since the body of 
the rule is false in /, U\ = {a}, U2 = {&}, and U3 = {a,b} are all nonempty 
unfounded sets in this case. It follows immediately by Theorem 12.91 that / 
is not a stable model. However, U2 is maximally /-consistent yet it is not a 
partial stable model because I* is not a minimal model of P 1 . 

Note that P has a unique (partial) stable model, which is ({&}, {a})- □ 

The characterizations for partial and total stable models in terms of un- 
founded sets provide a powerful tool for establishing relationships between 
stable models and partial stable models. 

3 UNFOLDING PARTIALITY 

In this section, we first show a translation for a disjunctive program into 
another disjunctive program. We then prove that the translation preserves 
the semantics of partial stable models. This result allows us to compute 
the partial stable models of a program by computing the stable models of 
the translated program. Finally we address the problem of query answering 
under the translation. 

3.1 TRANSLATION 

Let P be a disjunctive program. In the following, we describe a translation 
of P into another disjunctive program Tr(P) such that the stable models of 
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Tr(P) correspond to the partial stable models of P. 

Let us introduce a new atom a' for each a G Hb(P). An atom a' is 
said to be marked, and an ordinary atom a is then said to be unmarked. 
The intuitive reading of a* is that a is potentially true. For a set of literals 
L C Hb(P) U ~Hb(P), we define L* = {a' | a G L} U {~a' | ~a G L}. The 
translation Tr(P) of a disjunctive program P is as follows: 

Tr(P) = {A <- B, ~C"; A* <- B',~C\A^ B, ~C G P} U 

{a* <- a|a G Hb(P)} (2) 

where semicolons are used to separate program rules. Note that the Herbrand 
base of Hb(Tr(P)) is Hb(P) UHb(P)*. The rules a" <— a introduced for each 
a G Hb(P) enforce consistency in the sense that if a is true, then a must also 
be potentially true. 

Remark 3.1 Although for presentational purposes the translation is defined 
for ground programs, exactly the same translation applies to non-ground 
programs as well: for each predicate p we introduce a new predicate p*, 
hence for a (ground or non-ground) atom <f> = p{t\, . . . ,t„), the new atom 
is (p* = p*(ti, . . . , t n ) (cf. Example I3.5j) . Since our proofs do not depend 
on the assumption that a given program is finite, the conclusions reached 
cover also any non-ground program with function symbols whose semantics 
is determined by treating the program as a shorthand for its (possibly infinite) 
Herbrand instantation. 

A partial stable model of a given program will be interpreted by a corre- 
sponding stable model of the transformed program. The extra symbol a' for 
each atom a provides an opportunity to represent undefined (in three-valued 
logic) in terms of truth values of a* and a in two- valued logic. For each pair a 
and a', either of which can be true or false, there are four possibilities: when 
a' and a are in agreement, that is when they are both true or both false, 
the truth value of a is their commonly agreed truth value; the combination 
where a is false and a' is true then represents that a is undefined; and the 
fourth possibility where a is true and a* is false is ruled out by any models 
due to the consistency rules. This intended representation of a partial stable 
model is given by the following equations. 

Definition 3.2 Let M be a partial interpretation of a program P and N 
a total interpretation o/Tr(P). The interpretations M and N are said to 
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satisfy the correspondence equations if and only if the following equations 
hold. 



M t = {ae Hb(P) | a G N* and a' G N*} (CE1) 

M f = {ae Hb(P) | a G iV f and a' G iV f } (CE2) 

M u = {a G Hb(P) | a G iV f and a* G iV*} (CE3) 

= {a G Hb(P) laeiV* and a* G N f } (CE4) 

Note that total interpretations that are models of Tr(P) satisfy CE4 
immediately, since the set of rules {a* <— a | a G Hb(P)} is included in Tr(P). 
Consequently, the "fourth truth value" is ruled out. The following example 
demonstrates how the representation given in Definition 13.21 allows us to 
capture the partial stable models of a disjunctive program P with the total 
stable models of Tr(P). 

Example 3.3 Consider a disjunctive program 

P = {a V b <— ~c; b <— c <— ~c}. 

Now a becomes false by the minimization of partial models, since the falsity 
of a does not affect the satisfiability of any rule. Thus the unique par- 
tial stable model of P is M = (0, {a}). Note that the reduction Pm = 
{a V b <— u; b <— u; c <— u}. Then consider the translation 

Tr(P) ={aVb^ ~ c *; 
a' <— a; 

The unique stable model of Tr(P) is iV = {6*, c*} which represents (by CE2 
and CE3) the setting that b and c are undefined and a is false in M. □ 

It is well-known that a disjunctive program P may not have any partial 
stable models. In such cases, the translation Tr(P) should not have stable 
models either, if the translation Tr(P) is to be faithful. 

Example 3.4 Consider a disjunctive program 

P = {aV6Vc^; a <— b <— ~c; c <— ~a} 



b^ 


- 


c <— ~c*; 




— ~6; 


c* <— ~c; 


V * 


-6; 


c* <- c } 
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and its translation 

Tr(?) = {aV&Vc^; a< b'; b< c"; c< a*; 

a* V b* V c" <-; a* < 6; 6* <- ~c; c" < a } U C 

where C = {a* <— a; 6* 6; c* <— c} is the set of consistency rules. 

Consider a partial model M = ({a, 6},0) of P and a total model A" = 
{a,a',b,b*,c'} of Tr(P) that satisfy the equations CE1-CE4 in Definition 
13.21 Now the reduced program Pm is 

{aVfeVc^; a <— f ; 6 <— u; c^f} 

and since M' = ({a, b}, {c}) < M is a partial model of Pm, M is not a partial 
stable model of P. On the other hand, the reduct 

Tr{P) N = {aV6Vc^; a'V&'Vc'^ 6V}uC. 

But iV' = {a, a', 6, 6*} C iV is a model of Tr(P) N so is not a stable model 
of Tr(P). The reader may analyze the other candidates in a similar fashion. 
It turns out that P does not have partial stable models. Nor does Tr(P) 
have stable models. □ 

Partial stable models can be viewed as a logic programming account of 
the solution of semantic paradoxes due to Kripke In this account, un- 
defined means unknown for some individuals which will not lose semantic 
interpretations for other individuals. 

Example 3.5 Consider the following program with variables: 

P = { shave(bob, x) <— ~shave(a;, x)\ 

pay_by_cash(?/, x) V pay_by_credit(y, x) <— shave(x, y); 
accepted (x, y) <— pay_by_cash(x, y); 
accepted (x, y) <— pay_by_credit(x, y) }. 

This program intuitively says that Bob shaves those who do not shave them- 
selves; if x shaves y then y pays x by cash or by credit; either way is accepted. 
The predicate accepted is used here to demonstrate disjunctive reasoning. 

Assume there is another person, called Greg. Then clearly, we should 
conclude Bob shaves Greg, and Greg pays Bob by cash or by credit, either 
way is accepted. However, the program has no stable models in this case due 
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to the paradox whether Bob shaves himself or not. But it has two partial 
stable models, in both of which shave(greg, greg) is false and shave(bob, bob) 
is undefined (unknown). By translating the first two rules of P we obtain 

shaves(bob, x) <— ~shaves*(x, x); 
shaves* (bob, a;) <— ~shaves(x, x); 

pay_by_cash(?/, x) V pay_by_credit(y, x) <— shaves(a;, y); and 
pay_by_cash*(y, x) V pay_by_credit*(y, x) <— shaves* (x, y). 

The full translation Tr(P) yields a Herbrand instantiation over the universe 
{bob, greg} which has four total stable models. One of them is 

N = { shaves*(bob, bob), 

shaves(bob, greg), shaves*(bob, greg), 
pay_by_cash(greg, bob), pay_by_cash*(greg, bob), 
pay_by_credit*(bob, bob), accepted*(bob, bob), 
accepted (greg, bob), accepted* (greg, bob) }. 

Hence the fact that shaves(bob, bob) is undefined in the corresponding partial 
stable model M (recall the equations in Definition 13 .2j) is represented by 
shaves* (bob, bob) being true and shaves(bob, bob) being false in N. □ 

3.2 CORRECTNESS OF THE TRANSLATION 

The goal of this section is to establish a one-to-one correspondence between 
the partial stable models of a disjunctive program P and the (total) stable 
models of the translation Tr(P). It is first shown that the correspondence 
equations CE1-CE4 given in Definition 13 . 21 provide a syntactic way to trans- 
form a partial stable model M of P into a total stable model N of Tr(P) 
and back. More formally, we have the following theorem in mind. 

Theorem 3.6 Let M be a partial interpretation of a disjunctive program P 
and N a total interpretation of the translation Tr(P) such that CE1-CE4 are 
satisfied. Then M is a partial stable model of P if and only if N is a (total) 
stable model o/Tr(P). 

Our strategy to prove Theorem 13.61 is as follows: first, in two separate 
lemmas, we show the correspondence, in each direction, between unfounded 
sets for P and Tr(P) under M and N, respectively. These two lemmas are 
interesting in their own right as they show very tight conditions under which 
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the two previously studied notions of unfoundedness [HJ EEj are related. 
These results will then be used in the proof of the theorem. We first state 
two relatively simple facts. The first says that that the GL-transform has no 
effect on unfoundedness, and the second states that the translation preserves 
models via CE1-CE4 in Definition I3~2l 

Proposition 3.7 Let P be a disjunctive program and N a total interpreta- 
tion for P. Then, X C Hb(P) is an unfounded set for P w.r.t. N if and only 
if X is an unfounded set for P N w.r.t. N . 

PROOF. Note that A <- B G P N if and only if there is a rule A <- B, ~C G 
P such that C C X f , i.e., CniV t = |. Then it holds for any X C Hb(P) 
that 

X is not an unfounded set for P w.r.t. N 
3A^ P, ~C G P such that (1) Anl^t), (2) B n iV f = 0, 
(3) C n X 1 = 0, (4) P n X = 0, and (5) (A - X) n X 1 = 
3A «- P G P N such that (6) A n X ^ 0, (7) P n iV f = 0, 
(8) P n X = 0, and (9) (A - X) n X 1 = 
X is not an unfounded set for P N w.r.t. X. 

□ 




Proposition 3.8 Let M be a partial interpretation for a disjunctive program 
P and N a total interpretation for the translation Tr(P). Assume M and N 
satisfy the CEs. Then, M is a partial model of P if and only if N is a total 
model o/Tr(P). 

PROOF. It follows by the correspondence equations CE1-CE4 in Definition 
13.21 that M is not a partial model of P if and only if 

3 A <- P, ~C G P: M(V A) < M(B U ~C) 
<^ ~C G P: M (V A) < t and M(B U ~C) = t, or 

3 A «- P, ~C G P: M(\/ A) = f and M(P U ~C) = u 
<^ 3 A «- P, ~C" G Tr(P): X(V A) = f and X(P U ~C") = t, or 

3 A* <- P',~C G Tr(P): N(\J A') = f and X(P' U ~C) = t 
« 3 A <- P,~C" G Tr(P): X(\/A) < X(PU~C"), or 

3 A* <- P',~C G Tr(P): N(\J A') < X(P* U ~C) 

which is equivalent to stating that X is not a total model of Tr(P), since the 
consistency rules in Tr(P) are automatically satisfied by CE4. □ 
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Still assuming the setting determined by CEs, the following lemma gives a 
condition under which the unfounded sets w.r.t. N for Tr(P) can be converted 
into unfounded sets w.r.t. M for P. 

Lemma 3.9 Let P be a program, M a partial interpretation of P and N a 
total interpretation of the program Tr(P) such that CE1-CE4 are satisfied. 
Then, for any unfounded set X for Tr(P) w.r.t. N , the set of atoms Y = 
{a G Hb(P) | a* G X} is an unfounded set for P w.r.t. M. In addition, if X 
is N -consistent, then Y is M -consistent. 

PROOF. Consider any rule A <— B,~C G P such that A n Y ^ 0. It is 
proven in the sequel that one of the unfoundedness conditions UF1-UF3 ap- 
plies to A <— B, ~C. Two cases arise depending on the value of M(B U ~C). 

I. If M(B U ~C) = f , then UF1 is directly applicable. 

II. Suppose that M(B U ~C) ^ f which implies N(B' U ~C) = t by the 
CEs. Now A n Y ^ and the definition of Y imply A'nl^l. Since 
A' <— B*, ~C G Tr(P), X is an unfounded set for Tr(P) w.r.t. N and 
UF1 is not applicable to A* <— P*,~C, we know that either UF2 or 
UF3 applies to A' <- B*, ~C. 

(i) If UF2 applies to A' <- P', ~C, then P'nX ^ 0. It follows by the 
definition of Y that Bfl7^, i.e., UF2 applies to A <- P, ~C. 

(ii) If UF3 applies to A* <- P', ~C, then A^(V (A' - X)) = t. Since 
A* — X = (A — Y)' by the definition of Y, we obtain by the CEs 
that M(V (A - Y)) ^ f . Thus UF3 applies to A <- P, ~C. 

The proof of the consistency claim follows. To establish the contrapositive 
of the claim, suppose that Y is not M-consistent. Then Y fl M* 7^ 0, i.e., 
there exists an atom a G Hb(P) such that a E Y and a G M t . The former 
implies a' 6 I by the definition of Y while the latter gives us a 9 G N l by 
the CEs. Thus X n iV* 7^ and X is not iV-consistent. □ 

The next lemma shows that, under the specified conditions, an unfounded 
set for a given disjunctive program P corresponds to a collection of unfounded 
sets for the translation Tr(P). 

Lemma 3.10 Let M be a partial model of a disjunctive program P and N 
a total interpretation o/Tr(P) satisfying the CEs. If X is an M-consistent 
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unfounded set for P w.r.t. M, then Y = F U U where F = {a, a' \ a G X} 
and U C {a | a G X f , a* G X*} is an unfounded set for Tr(P) w.rJ. X. 

PROOF. Let X be an M-consistent unfounded set for P w.r.t. M and let 
F = F U P satisfy the requirements above. Since any atom in F is either 
marked or unmarked, two cases arise. 

I. Suppose that a* G Y which implies by the definition of Y that a 6 7 . 
Then it is clear that that UF2 applies to the consistency rule a' <— 
a G Tr(P). Let us then prove that one of the unfoundedness conditions 
applies to any rule A' <— B', ~C G Tr(P) satisfying a' G A'. Since 
X is a total interpretation, we have N(B* U ~C) = f (in which case 
UF1 applies to A' *- B*,~C) or N(B m U ~C) = t in which case 
M (B U ~C) > f . Since a G X and a G A, and UF1 does not apply to 
A <— P, ~C, we only need to consider UF2 and UF3. If UF2 applies 
to A <— B, ~C, 3b G P such that 6 G X. It follows by the definition of 
Y that 6* G F. Hence UF2 applies to A' «- P*, ~C If UF3 applies to 
A <- P, ~C, 36 G A such that M(b) > f and b £" X. Then we know 
that X(6') = t by the CEs. Further, by the definition of Y, b g X 
implies 6* g F. Hence UF3 applies to A* «- P', ~C. 

II. Suppose that a G F. Then consider any rule A P,~C* G Tr(P) 
such that a G A. Since X is a total interpretation, N(B U ~C*) = f 
(in which case UF1 applies to A <- P, ~C") or X(P U ~C*) = t. In 
the latter case, we know that 36 G A such that N(b) = t, since X is a 
model of Tr(P) by Proposition 13.81 (recall that M is a partial model of 
P). Then suppose that b G F, i.e., 6 G P or 6 G U by the definition of 
F. If 6 G P, then 6 G X by the definition of P. On the other hand, 
N(b) = t implies M(b) = t. Thus M t n X ^ 0, contradicting the 
M-consistency of X. If 6 G P, the the definition of U implies N(b) = f, 
a contradiction. Hence, b G" F and UF3 applies to A <— P, ~C*. 

□ 

We note that the M-consistency of X is also a necessary condition for 
the correspondence to hold. 

Example 3.11 Consider a disjunctive program P = {aV6^; a <— ~a} 
and a partial model M = ({b}, 0) of P. It can be checked easily that X = {b} 
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is an unfounded set for P w.r.t. M: UF3 applies to the only rule in which b 
appears in the head. But X is not M-consistent. Now consider 

Tr(P)={aV&^; a'V6*^; a< a'; 

a* <— ~a; a' <— a; b' <— b }. 

The total interpretation corresponding to M above is N = {a',b,b'}. How- 
ever, Y = {a, b, b'} is not unfounded for Tr(P) w.r.t. N, since for b G Y and 
the first rule in Tr(P), none of the unfoundedness conditions applies. □ 

Let us establish Theorem 13.61 in two separate theorems. 

Theorem 3.12 Let P be a disjunctive program. If N is a stable model of 
the translation Tr(P), then the partial interpretation M of P satisfying the 
correspondence equations CE1-CE4 is a partial stable model of P. 

PROOF. Let N be a stable model of Tr(P). Then it follows by the pres- 
ence of consistency rules {a* <— a | a G Hb(P)} in Tr(P) that there is no 
a G Hb(P) such that a G iV and a' G" N, since iV is a total model of 
Tr (P) . Thus it makes sense to define M as the partial interpretation satisfy- 
ing CE1-CE4. We prove that M t is a minimal total model of P M , and M f 
is a maximal M-consistent unfounded set for P w.r.t. M. 

I. Let us first establish that for any rule A <— B, ~C G P, A <— B G P M 
A «- B G Tr(P) N . So consider any A <- P, ~C G P. It 
follows by the CEs and the definitions of P M , Tr(P) and Tr(P) N that 
Af-Be P M there is a rule A «- P, ~P G P such that P C M f 

there is a rule A «- P, ~P* G Tr(P) such that D* C A^ f 
A <— B E Ti{P) N . Note that within these equivalences A <— P, ~C 
and A <— P, ~P need not be the same rules of P. 

II. Let us then prove that M t is a minimal total model of P M . If we 
assume the contrary, two cases arise. 

- M* is not a total model of P M , i.e., there is a rule A <— B G P M 
such that ilP(P) = t, but M*(A) = f. It follows by the CEs that 
N(B) = t and N(A) = f . Thus A «- P is not satisfied in and 
thus iV is not a model of Tr (P) N , as P M C Tr (P) N holds by (I) 
above. A contradiction, since N is a stable model of Tr(P). 
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- There is a total model M' of P M such that M' C M t . Then define 
a total interpretation AT' = M' U {a* | a* G N}. ByM'c M* and 
the CEs, we obtain N f C N (only some unmarked atoms of N 
are not in N'). Since M' is a total model of P M , and M' and A 7 
coincide on the atoms of Hb(P), every rule in P M is satisfied by 
N'. By (I), the difference Tr(P) N — P M contains only consistency 
rules a' <— a (for every a G Hb(P)) and rules of the form A* <— P* 
(for some A <— P, ~C G P). These rules are all satisfied by N', 
since N is a total model of Tr(P) N , N' C N, and AT' and AT 
coincide on the marked atoms in Hb(P)\ Thus N' is a total 
model of Tr(P) . Then N' C N implies that iV is not a minimal 
model of Tr(P) N nor a total stable model of P. A contradiction. 

III. Since N is a total stable model of Tr(P), it holds by Theorem 12.81 
that A^ f is the greatest unfounded set for Tr(P) w.r.t. N. More- 
over, N f is N-consistent, since iV'niV' = 0. Note that a* G A^ f 
implies a G A^ f , since N satisfies a* <— a G Tr(P). Thus M f = 
{a G Hb(P) \aeN f and a' G A^ f } = {a G Hb(P) | a' G A^ f }. It fol- 
lows by Lemma 13.91 that M f is an M-consistent unfounded set for P 
w.r.t. M. 

Then assume that M f is not maximal, i.e., there is an M-consistent 
unfounded set X for P w.r.t. M such that X D M f . So there is an 
atom a G X such that a ^ M f . Then a ^ M f implies a G M t or 
a G M u . In both cases, by the CEs, a' G AT*, i.e., a' ^ A^ f . Then 
construct F = {a, a' | a 6 X}. According to Lemma UllOl that X is 
an M-consistent unfounded set for P w.r.t. M implies that F is an 
unfounded set for Tr(P) w.r.t N. However, a & X implies a' G Y 
but a* ^ N f . Thus a* G AT* indicating that X is not unfounded free 
for Tr(P). Consequently, N is not a stable model of Tr(P) by the 
characterization of stable models in Theorem 12.81 a contradiction. 

□ 

Theorem 3.13 Let P be a disjunctive program. If M is a partial stable 
model of P , then the total interpretation N satisfying the correspondence 
equations CEl-CEJf. is a stable model of the translation Tr(P). 

PROOF. Suppose that M is a partial stable model of P. Then we know by 
Theorem 12.91 that (i) M t is a minimal total model of P M and (ii) M f is a 
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maximal M-consistent unfounded set for P w.r.t. M. Then define N as the 
total interpretation of Tr(P) satisfying the CEs. It follows by Lemma fc. 101 
that X f = {a,a'\ae M f } U M u is an unfounded set for Tr(P) w.r.t. N. 

Let us then assume that N is not a stable model of Tr(P). Equivalently, 
it holds by Theorem 12.81 that X f is not the greatest unfounded set for Tr(P) 
w.r.t. N. So there is an unfounded set X for Tr(P) w.r.t. X such that 
X f C X and X* n X ^ hold. It follows by Proposition EHI that X is also 
an unfounded set for Tr(P) N w.r.t. X. 

Then consider any A <— P G P M for which there is a rule A <— P, ~C G P 
such that C C M f . It follows by the CEs that C C N f . Since A «- P, ~C" G 
Tr(P), it follows that A <- P G ^(P)^. Thus P M C ^(P)^ holds, as 
Tr(P) N contains among others the consistency rules {a* <— a \ a G Hb(P)}. 

Recall that M t = JV* n Hb(P) is a minimal total model of P M . We also 
distinguish a set of atoms X' = X n Hb(P). Let us then establish that X' is 
an unfounded set for P M with respect to M* in the two-valued sense. 

I. If X' is not such a set, it follows by Definition l2.3l that there is A <— B G 
P M with Af]X' ^ such that BCM\ Bf]X' = and (A-X')nM t = 
0. It follows that A ^ B e Tr(P) N , as P M C Ti(P) N . Since A and 
P are subsets of Hb(P), M t = AT* n Hb(P) and I' = ln Hb(P), we 
obtain A n X ^ 0, P C N\ B n X = and (A - X) n iV* = 0. Then 
X is not an unfounded set for Tr(P) N w.r.t. X, a contradiction. 

It follows by Lemma El that M* - X' is a total model of P M . It follows 
by the minimality of M t that M t n X' = and M* n X = 0. Moreover, it 
follows by Lemma EHl that F = {a G Hb(P) | a* G X} is an unfounded set 
for P w.r.t. M. It remains to establish that Y is M-consistent and M f C Y. 

II. Suppose that Y is not M-consistent, i.e., it holds for some a G Hb(P) 
that (a) a G F and (b) a G M*. Now (b) implies by the CEs that 
a G X* and a* G X*. On the other hand, it follows by (a) and the 
definition of Y that a* G X. Thus one of the unfoundedness conditions 
applies to the rule a' <— a G Tr(P), as X is an unfounded set for 
Tr(P) w.r.t. X. Now UF1 is not applicable, as a ^ X f , and UF3 is not 
applicable, as a* G X. Thus UF2 must be applicable to a' <— a. It 
follows that a G X, too. Then there is a G Hb(P) such that a G M t 
and a G X. A contradiction with M* PI X = established above. 

III. Consider any a G M f . Thus a G Hb(P) and a' G X f follows by the 
CEs. Then X f C X implies a' G X as well as a G Y by the definition 
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of F. Thus M f C Y. On the other hand, recall that JV* n X ^ and 
M t n X = 0. Then a* G iV* n X holds for some a G Hb(P). It follows 
that a' G iV* and a* G X. The former implies a ^ M f by the CEs. 
The latter implies a G F by the definition of F. Hence M f C F. 

Thus M f is not a maximal M-consistent unfounded set for P w.r.t. M, a 
contradiction. Hence N must be a stable model of Tr(P). □ 

It is worthwhile at this point to briefly comment on the proof of The- 
orem ESI as given in (42], which proceeds in several steps. Given two par- 
tial interpretations M and M 1 of a disjunctive program P, let N and N' , 
respectively, be the corresponding total interpretations of Tr(P) such that 
CE1-CE4 are satisfied. Firstly, it can be shown that M is a partial model 
of Pm> if and only if N is a total model of Tr(P) . Secondly, since the 
truth-ordering for partial interpretations corresponds to the subset ordering 
for total interpretations, it can be shown that the minimal partial models of 
Pm> correspond to the minimal total models of Tr(P) N . Thirdly, based on 
a characterization of partial models in general jlH], we conclude that M is a 
partial stable model of P if and only if N is a total stable model of Tr(P). 

Looking back to results established so far, we know by Theorem 13. 131 that 
any partial stable model M of P can be mapped to a stable model 

f(M) = M* U (M* U M u ) # (3) 

of Tr(P). Similarly, any stable model N of Tr(P) can be projected to a 
partial stable model 

g(N) = ({a G Hb(P) | a G N}, {a G Hb(P) | a N and a* £ N}} (4) 

of P by Theorem 13.121 These equations and the corresponding theorems 
indicate that / and g are functions between the set of partial stable models of 
P and the set of stable models of Tr(P). In the sequel, it is established that 
these functions are bijections, which means that our translation technique 
does not yield spurious models for programs although new atoms are used. 
This is highly desirable from the knowledge representation perspective. 

Theorem 3.14 The partial stable models of a disjunctive program P and 
the total stable models of the translation Tr(P) are in a one-to-one corre- 
spondence. 
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PROOF. Let / and g be defined by the equations © and PJ, respectively. 
It is straightforward to see that / is injective, i.e., /(Mi) = f(M 2 ) implies 
Mi = M 2 . Then assume that g(Ni) = g(N 2 ) holds for some stable models Aq 
and N 2 of Tr(P). It follows by the definition of g for any a G Hb(P) that (i) 
a G iVi a e N 2 and (ii) a G" Aq and a* G" N\ a N 2 and a' N 2 . 

Then consider any a G Hb(P) such that a* G i\q. Two cases arise. If a G Aq, 
it follows by (i) that a G A 2 . Since N 2 satisfies the rule a' <— a G T^P)^ 2 , 
we obtain a* G A 2 . On the other hand, if a G" Aq it follows by (i) that a G" A 2 . 
Assuming that a* G" N 2 implies by (ii) that a' G" Aq, a contradiction. Hence 
a* G iV 2 also in this case. By symmetry, a' G N 2 implies a' G Aq. 

Thus it holds for any a G Hb(P) that (iii) a' G JVi a' G iV 2 . It 

follows by (i) and (iii) that Ni = N 2 so that g is injective, too. Thus / and g 
are bijections and inverses of each other, as g(f(M)) = M and f(g(N)) = N 
hold for any (partial) stable models M and N. Hence the claim. □ 

None of the preceding proofs relies on the assumption that the given 
program is finite. Therefore, all of these results presented in this section 
apply in the non-ground case as well. 

3.3 QUERY ANSWERING 

Let us yet address the possibility of using an inference engine for comput- 
ing total stable models to answer queries concerning partial stable models. 
This is highly interesting, because there are already systems available for 
computing total stable models [U [211 EH HE] while partial stable models lack 
implementations. Here we must remind the reader that partial stable models 
can be used in different ways in order to evaluate queries. Typically two 
modes of reasoning are used: certainty inference and possibility inference. 
In the former approach, a query Q should be true in all (intended) models 
of P while Q should be true in some (intended) model of P in the latter 
approach. Moreover, maximal partial stable models (under set inclusion) are 
sometimes distinguished; this is how regular models and preferred extensions 
are obtained for normal programs [101123 EH) • We are particularly interested 
in possibility inference where the maximality condition makes no difference 
(see [iniQHl for certainty inference): M(Q) = t for some partial stable model 
M of P if and only if M'(Q) = t for some maximal partial stable model M' 
of P. 

We consider queries Q that are sets of literals over Hb(P) and queries 
are translated in harmony with the CEs: Tr(Q) = Q U Q* . As a direct 
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consequence of Theorem 13. 131 and CE1, we obtain the following. 



Corollary 3.15 A query Q is true in a (maximal) partial stable model of P 
if and only ifTr(Q) is true in a stable model ofTr(P). 

What about using a query answering procedure for partial stable models 
to answer queries concerning stable models? A slight extension of the trans- 
lation Tr(P) is needed for this purpose: let Tr 2 (P) be Tr(P) augmented with 
a set of rules {/ <— a*, ~a | a e Hb(P)} where / ^ Hb(P) is a new atom. The 
purpose of these additional rules is to detect partial stable models with re- 
maining undefined atoms. A query Q is translated into Tr 2 (Q) = Q U {~/}. 

Corollary 3.16 A query Q is true in a stable model of P if and only if 
Tr 2 (Q) is true in a partial stable model o/Tr 2 (P). 

This result allows query answering for stable models to be conducted by 
a procedure for partial stable models, e.g., by the abductive procedure of 
Eshghi and Kowalski [T6] . 

4 UNFOLDING DISJUNCTIONS 

In this section we develop a method for reducing the task of computing a 
(total) stable model of a disjunctive program to computing stable models for 
normal (disjunction-free) programs. This objective demands us to unfold 2 
disjunctions from programs in a way or another. Since the problem of de- 
ciding whether a disjunctive program has a stable model is S 2 -complete [T2] 
whereas the problem is NP-complete in the non-disjunctive case [Hi], the 
reduction cannot be computable in polynomial time unless the polynomial 
hierarchy collapses. This is why our reduction is based on a generate and 
test approach. 

The basic idea is that given a disjunctive program P we compute its 
stable models in two phases: (i) we generate model candidates and (ii) test 
candidates for stability until we find a suitable model. For generating model 

2 The idea of unfolding disjunction generally refers to performing some transformations 
on disjunctions in order to remove them ^flElE]. However, such transformations do not 
necessarily remove all disjunctions or do not preserve stable semantics. 
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candidates we construct a normal program Gen(P) such that the stable mod- 
els of Gen(P) give the candidate models. For testing a candidate model M 
we build another normal program Test (P,M) such that Test (P,M) has no 
stable models if and only if M is a stable model of the original disjunctive 
program P. Hence, given a procedure for computing stable models for nor- 
mal programs all stable models of a disjunctive program P can be generated 
as follows: for each stable model M of Gen(P), decide whether Test (P,M) 
has a stable model and if this is not the case, output M as a stable model of 
P. This kind of a generate and test approach is used also in dlv [21] which 
is a state-of-the-art system for disjunctive programs. The difference is that 
we reduce the test and generate subtasks directly to problems of computing 
stable models of normal programs whereas in dlv special techniques for the 
two subtasks have been developed based on the notion of unfounded sets for 
disjunctive programs. 

It is easy to construct a normal program for generating candidate models 
for a disjunctive program P. Consider, e.g., a program G0(P) which contains 
for each atom a G Hb(P), two rules a <— ~d; a <— ~a where d is a new 
atom denoting the complement of the atom a, i.e., a is in a stable model 
exactly when a is not. These rules generate stable models corresponding to 
every subset of Hb(P). In order to prune this set of models to those with all 
rules in P satisfied, it is sufficient to include a rule 

/ <- ~/, . . . , ~a fc , bx,...,b m , ~c x , . . . , ~c n (5) 

for each rule of the form (QJ) in P where / is a new atom. As / cannot be in 
any stable model, the rule functions as an integrity constraint eliminating the 
models where each foj is included, every Cj is excluded but no a\ is included. 

In order to guarantee completeness, it is sufficient that for each stable 
model M of P there is a corresponding model candidate which agrees with 
M w.r.t. Hb(P). It is clear that G0(P) satisfies this condition. However, for 
efficiency it is important to devise a generating program that has as few as 
possible (candidate) stable models provided that completeness is not lost. An 
obvious shortcoming of G0(P) is that it generates many candidates even if 
the program P is disjunction-free. In order to solve this problem we construct 
for given a disjunctive program P a generating program G1(P) as follows: 

G1(P) = {a< a,B,~C | A *- B,~C G P D ,a G A} U 

{d <— ~a | a G Heads(Pn)} U 

{/ < / ; „a, B, ~C \ A<— B, r^c G P D } U P N 
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where Pn is the set of the normal rules in P and Pd are the other (proper 
disjunctive) rules in P, i.e. P = Pn U Pd, and Heads (Pd) is the set of atoms 
appearing in the heads of the rules in Pd. 

The program G1(P) has typically far fewer stable models than G0(P) and 
the number of "extra" candidate models which do not match stable models 
of P is related to the number of disjunctions in P. For example, if P is 
disjunction-free, the stable models of G1(P) correspond exactly to the stable 
models of P. However, for a disjunctive program P, G1(P) can easily have 
"extra" stable models. Consider, e.g., 

P = {aV6^} (6) 

for which G1(P) = {a <— ~a; a <— ~a; b <— ~6; b <— / <— 

~/, ~a, ~6} has a stable model, {a, 6}, not corresponding to a stable model 
of P. In fact, G1(P) only requires for each proper disjunctive rule in P that 
some non-empty subset of the head atoms of the rule is included in the model 
candidate when the body of the rule holds. Hence, for such a rule with d 
disjuncts in the head there are 2 d — 1 possible subsets and in the worst case 
2 d — 2 of these could lead to "extra" model candidates. This means that 
in the worst case G1(P) can have an exponential number of "extra" model 
candidates w.r.t. the number of disjunctions in P. 

In order to decrease the number of "extra" models we introduce a tech- 
nique exploiting a key property of supported models [3j: each atom a true 
in a model M of P must have a rule supporting it, i.e., there is a rule A <— 
B,~C E P such that a e A, M(B U ~C) = t, and M(V (A - {a})) = f . 
Since every stable model of P is also a supported model of P, it makes per- 
fect sense to require supportedness from the candidate stable models. For 
this, we introduce a new atom a s , which denotes the fact that atom a has 
a supporting rule, for each atom a appearing in the head of a disjunctive 
rule. The intuition behind the set of rules Supp(P) below is that a rule can 
support exactly one of its head atoms and we may exclude every model that 
has an atom without a supporting rule: 

Supp(P) = 

{a s < (A - {a}), B, ~C \ A^ B,~C e P,ae An Heads(P D )} U 

{/ «- ~/ 3 a > ~ flS I a ^ Heads(P D )} (7) 

where Heads (Pd) is the set of atoms appearing in the heads of the proper 
disjunctive rules in P. For example, for P in (JHJ), 
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Supp(P) = {a s <- 6 s <- ~a; / < /, a, ~a s ; / <- ~/, 6, ~6 S }. 

Now G1(P) U Supp(P) has exactly two stable models {a, a s , b} and {b, b s , a} 
corresponding to the two stable models {a} and {b} of P. 

Combining this idea with G1(P) gives a promising generating program 

Gen(P) = G1(P) U Supp(P) (8) 

which still preserves completeness. 

Proposition 4.1 Let P be a disjunctive program. Then if M is a stable 
model of P, there is a stable model N of Gen(P) = G1(P) U Supp(P) with 
M = iVnHb(P). 

PROOF. Let M be a stable model of P and 

N = M U {a | a G Heads(P D ) - M} U {a s | a G M n Heads(P D )}. 

Now clearly M = iVnHb(P). We show first that (i) iV is a model of Gen(P) 7V 
and then that (ii) if there is a model N' of Gen(Pj such that N' C N then 
AfCJV' holds. These together imply that N is a stable model of Gen(P). 

For property (i) consider rules in Gen(P) Ar = G1(P) N U Supp(P) JV start- 
ing with those in G1(P) N . Suppose a <— B e G1(P) N . If a G Heads(P D ), 
then a £ N and, hence, a G M C N. Otherwise if a G" Heads (Pd), then 
a <- P G P M implying that iVQ M) satisfies a <— B. If a «— G G^P)^, 
then a G" M and a & N. If / <- P G G^P)^, then there is a rule 
A <- P, ~C G P D such that AnM = and A ^ 5 G P M . As M is a model 
of P M , P ^ M and, consequently, P £ TV. Thus, iV is a model of Gl(P) N . 
Next consider rules in Supp(P) 7V . If a s <— P G Supp(P) Ar , then there is a 
rule A <- P, ~C G P such that A <- P G P M and (A - {a}) n M = 0. If 
BCiV, then PCM and, hence, AnM^0asMisa model of P M . Thus, 
a G M. If / <- a G Supp(P) A ', then a s G" TV, and a ^ N. This implies that 
TV is a model of Supp(P) and that (i) holds. 

For property (ii) consider a model N' of Gen(P) 7V such that iV' C iV. First 
we show that N' n Hb(P) is a model of P M implying that AfCiV'n Hb(P). 
Consider A <- B G P M . If the body P is true in TV'nHb(P) C iVnHb(P) = 
M, then at least one a G A n M. Then a <- P G Gen(P) Ar and a G iV'. 
Hence, MCiV'n Hb(P). UaeN, then a G" iV and, hence, a ^G G1(P) W 
implying a G iV'. If a s G iV, then a G M fl Heads(PD). Then there is a rule 
A <- B G P M such that PCM but (A — {a}) fl M = 0. This is because 
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otherwise M — {a} would be a model of P M contradicting the minimality of 
M. Hence, a s «- B G Supp(P) iV and B C M C JV' implying a s G iV'. Thus, 
N C N' and (ii) holds. □ 
A (total) model candidate M C Hb(P) is a stable model of a program 
P if it is a minimal model of the GL-transform P M of the program. This 
test can be reduced to an unsatisfiability problem in propositional logic using 
techniques presented in |3B]: M is a minimal model of P M if and only if 



P M U {->a | a G Hb(P) - M} U {-.&x V • ■ ■ V ^b m } (9) 

is unsatisfiable where M = {bi, . . . ,b m } and the rules in P M are seen as 
clauses. This can be determined by testing non-existence of stable models for 
a normal program Test(P, M) which is constructed for a disjunctive program 
P and a total interpretation M C Hb(P) as follows: 

Test(P, M) — {a <— ~a, B \ A ^ B e P^ 1 , a e An M,B <Z M} U 
{a <— ~a | a G Heads(PD)} U 

{/ < /, ~A, B | A <- B G B C M} U 

{a <- P G P^ 7 | a G M, P C M} U 

where Pn is the set of the normal rules in P and Pq are the proper disjunctive 
rules in P and Heads(Po) is the set of atoms appearing in the heads of the 
rules in Pd. The idea is that stable models of Test(P, M) capture models of 
the reduct P M that are properly included in M. 

Proposition 4.2 Let P be a disjunctive program and M a (total) model of 
P. Then M is a minimal model of P M if and only if Test (P, M) has no 
stable model. 

PROOF. Let M C Hb(P) be a total model of P. 

(=►) Let iV be a stable model of Test (P, M). If a G Hb(P)-M, then there 
is no rule with a in the head in Test(P, M) and a G" N '. Hence, N H Hb(P) C 
M. As / ^ and / <- M G Test(P, M)^, there is some a G M such that 
a ^ n Hb(P). Consider A ■<— B E P M . Let P C n Hb(P) C M but 
suppose An A^nHb(P) =0. If A = {a}, a <- P G Test(P,M) JV and a G A r , a 
contradiction. Otherwise / P G Test(P, M)^ and /6iV,a contradiction. 
Hence, AT n Hb(P) is a model of P M but A^ n Hb(P) C M implying that M 
is not a minimal model of P M . 
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(<=) Assume that M is not a minimal model of P M . As M is a model 
of P M , there is a minimal model M' C M of P M . We show that A = 
M' U {a | a G Heads(P D ) - M'} is a minimal model of Test(P, M) N , i.e., A a 
stable model of Test(P, M). Now Test(P, M) N = 

{a<— B\A<— Be Pjf, a <E Af) M f) M', B C M} U 

{a <— | a G Heads(P D ) - M'} U 

{/ <- B | A <- P G P^, inM' = |,BCM}U 

{fl^BeP N M |flGM,BCM}U 

{/ <— M}. 

It is easy to check that (i) A is a model of Test(P, M) N . Assume there is 
a model A' of Test(P,M) 7V such that A' C A holds. We show that then 
(ii) A C A' holds as follows. We notice that for all a G Hb(P), a G A' 
implies a G M' . Consider A <- P G P M . If P is true in A' n Hb(P), then P 
is true in M' and, thus, PCM and some aeinMnM'. Then a <— P G 
Test(P, M)^ and a G A'. Hence, A'nHb(P) is a model of P M which implies 
M' C A' n Hb(P). If a G A, then a £ M' and a ^G Test(P, M) N implying 
a G A'. Then A C A' holds. Now (i) and (ii) imply that A is a minimal 
model of Test(P, M) N and, hence, a stable model of Test(P, M). □ 

Example 4.3 Consider a disjunctive program P and its generator Gen(P): 

P = { a yb< c} 

Gen(P) = { a 
a 
f 



For a stable model {b, b s , a} of Gen(P) the corresponding model candidate is 
Mi = {b, b s , a} fl Hb(P) = {b} and the test program: 

Test (P,Mi) = {b < 6; 

a <— ~a; b <— ~6; 

/ <- ~/, ~a, ~6; / <- ~/, b } 

Test(P, Mi) has no stable models and, hence, Mi is a stable model of P. □ 



~/, ~a, ~6, ~c; 
~6, ~c; 6 s <— ~a, ~c; 
~/, a, ~a s ; / < /, b, ~6 S } 
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The simple generate and test paradigm can be optimized by building model 
candidates gradually. This means that we start from the empty partial in- 
terpretation and extend the interpretation step by step. An interesting ob- 
servation is that the technique for testing minimality can be used to rule out 
a partial model candidate of Gen(P) at any stage of the search and not just 
when a total model of the program P has been found. This can be done by 
treating a partial interpretation M as a total interpretation where undefined 
atoms are taken to be false and using the Test(P, M) program. 

Proposition 4.4 Let P be a disjunctive program and M a total interpreta- 
tion. // Test (P,M) has a stable model, then there is no (total) stable model 
M' ofP such thatM C M' . 

PROOF. Let Test(P, M) have a stable model. As shown in the proof of 
Proposition S3 then there is a model M" of P M with M" C M. Consider 
any total interpretation M' such that M C M' and M' is a model of P M . 
Now M' is not a minimal model of P M as P M C P M and, hence, M" is a 
model of P M ' but M" C M C M' . □ 

Notice that for a total interpretation M, Proposition 14.41 can only be 
used for eliminating stable models of P extending M. For guaranteeing the 
existence of a stable model of P, a total model of P needs to be found making 
Proposition 14.21 applicable. 

Our approach to testing minimality of model candidates differs from that 
used in DLV \22\ . We check minimality by directly searching for models of 
the reduct strictly contained in the candidate model. In DLV a dual approach 
is used based on the notion of unfounded sets for disjunctive programs [2H] 
and minimality testing is done using a SAT solver. Our approach could be 
implemented straightforwardly using a SAT solver, too, but we have chosen 
to use the same logic program core engine for generating and testing sub- 
tasks in order to keep the implementation as simple as possible. A basic 
difference is that in our approach the set of clauses (jUJ used for minimal- 
ity testing follow the structure of the original program whereas in the DLV 
approach dual clauses (with each literal complemented) are employed. More- 
over, DLV employs a couple of optimizations which have not been exploited 
in our approach. First, DLV adopts specialized algorithms for some syntac- 
tically recognizable classes of rules like head cycle free programs. Second, 
DLV employs modular evaluation techniques for minimality testing where the 
program is divided into components based on its dependency graph and the 
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minimality of a candidate model is tested for each component separately 
by exploiting specialized algorithms for components with corresponding re- 
stricted form whenever possible. For a more detailed comparison, see [22]. 

5 EXPERIMENTS 

In this section, we compare dlv [21], a state-of-the-art implementation of the 
stable model semantics for disjunctive logic programs, with an implementa- 
tion of the generate and test approach of the previous section which we call 
GnT. In Section I57T1 we explain briefly implementation techniques employed 
in GnT and explain the setup for the experiments. For comparisons we 
use three families of test problems related to reasoning about minimal mod- 
els ^2], evaluating quantified Boolean formulas [47J, and planning [34J for 
which encodings of the problem instances as logic programs and test results 
are presented in Sections f5.2H5.4| respectively. All benchmarks used in the ex- 
periments are available at http : //www . tcs .hut . f i/Sof tware/gnt /benchmarks/ jnssy- tests -20 

5.1 IMPLEMENTATION 

The implementation of GnT [H| is based on SMODELS [46, |41|, a program 
that computes stable models of normal logic programs. The basic idea behind 
GnT is to use two instances of the smodels engine, one that generates the 
model candidates and one that checks if they are minimal. To implement the 
idea it is enough to extend the SMODELS engine only slightly. Figure [T] shows 
the pseudo-code for GnT modified from the original smodels function pre- 
sented in [IH]. The function gnt(G, P, A) takes as input a normal (generator) 
program G, a disjunctive program P and a partial model (a set of literals) 
A and performs a backtracking search for stable models of G. It returns 
a stable model M of G which agrees with A and for which minimal(P, M) 
returns true if such a stable model exists and otherwise it returns false. It 
uses functions expand(G, A), extend(G, A), conflict(G, A), heuristic(G, A), 
and minimal(P, A). The first four are as in the original smodels procedure: 

• expand(G, A) returns a partial model which expands the given partial 
model A by literals satisfied by all (total) stable models of G agreeing 
with A (obtained using a generalized well-founded computation); 
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• extend(G, A) returns a partial model extending the partial model A by 
literals obtained by expand enhanced with lookahead techniques. 

• conflict(G, A) checks whether there is an immediate conflict, i.e., if the 
partial model A contains a complementary pair of literals and 

• heuristic(G, A) returns an atom undefined in A to be used as the next 
choice point in the backtracking search for stable models. 

For further details on these functions see [Hf. The function minimal (P, A) 
performs the minimality test for a disjunctive program P and a partial model 
A given in Proposition 14.41 using a call to SMODELS, i.e., it views A as a total 
model A' where all atoms undefined in A are taken to be false, builds the 
program Test(P, A'), calls SMODELS and returns false if Test(P, A') has a 
stable model and otherwise returns true. To compute a stable model for a 
disjunctive program P, the procedure ^nt(Gen(P), P, 0) is called. First gnt 
extends the given partial model and checks for conflicts. If all atoms are 
covered by the extended partial model, then a (total) model candidate has 
been found and it is checked for minimality. Otherwise the heuristic function 
selects a new undefined atom x and gnt searches recursively first for models 
where x is false. If no such model is found, the partial model is expanded 
by making x true. If there is a conflict or the expanded model does not 
pass an "early" minimality test, the procedure backtracks and otherwise it 
continues the search recursively using the expanded model. As the "early" 
minimality tests are computationally quite expensive, some optimization has 
been employed so that such tests are performed only when backtracking from 
a model candidate. For this there is a global variable 'WasCovered' which 
is initially set to false and which is set to true when a model candidate 
is found. However, it should be noticed that when backtracking from a 
model candidate, the test could be repeated at each backtracking level until 
it succeeds. The implementation of the gnt procedure shown in Figure Q] 
consists of a few hundred lines of code |H] on top of the SMODELS system. 

In the sequel, we report several experiments which we carry out in order 
to compare dlv (version 2003-05-16) with GnT which is based on SMOD- 
ELS (version 2.27) and uses lparse (version 1.0.13) as an instantiator. We 
consider two versions of our approach, GnTI and GnT2, which are simi- 
lar except that in GnTI generating program G1(P) is used and in GnT2, 
Gen(P) = G1(P) U Supp(P). All of our tests are run under Linux 2.4.20 
operating system on a 1.7 GHz AMD Athlon XP 2000+ computer with 1 
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function gnt(G, P, A) 
A := extend(G,A) 
if conflict(G, A) then 

return false 
else if A covers Hb(G) then 

WasCovered := true 

if minimal(P, A) then 
return A 

else 

return false 

end if 
else 

x := heuristic(G, A) 
A' := gnt(G, P,AU 
if A' 7^ false then 

return A' 
else 

A' := expand(G, A U {x}) 
if conflict (G, A') then 

return false 
else if WasCovered then 

if not minimal(P, A') then 
return false 

end if 
end if 

WasCovered := false 
return gnt(G,P,A') 
end if 
end if. 

Figure 1: GnT Procedure 
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GB of memory. Execution times are measured using the customary Unix 
/usr/bin/time command. 

5.2 MINIMAL MODELS 

Our first test problem is the Sf-complete problem of deciding the existence of 
a minimal model of a set of clauses in which some specified atoms are true [12J. 
This problem is mapped to a stable model computation problem as follows. 
For a problem instance consisting of a set of clauses and some specified atoms, 
a program P is constructed where each clause a\ V • • • V a n V -161 V • • • V -b m 
is translated into a rule a% V - • • Va n <— bi, . . . , b m and for each specified atom 
Cj, a rule 

/<-~/,~c< (10) 

is included. Now P has a stable model if and only if there is a minimal model 
of the clauses containing all specified atoms q. 

The test cases (random disjunctive 3-SAT programs) are based on random 
3-SAT problems having a fixed clauses/atoms ratio c and they are constructed 
as follows. Given a number of atoms n, a random 3-SAT problem is gener- 
ated, i.e. c x n clauses are generated each by picking randomly three distinct 
atoms from the n available and selecting their polarity uniformly. This is done 
using a program makewff developed by Bart Selman. Then the clauses are 
translated into rules as described above and for % = 1, . . . , |_2n/100j and for 
random atoms q, the extra rules ([TO are added. The problem size is con- 
trolled by the number of atoms n which is increased by increments of 10. 
For each n, we test 100 random 3-SAT programs and measure the maximum, 
average, and minimum time it takes to decide whether a stable model exists. 

In the first set of tests we study the effect of different generating programs 
on the performance of our approach, i.e., we compare GnTI and GnT2, 
which are similar except that in GnTI generating program G1(P) is used 
and in GnT2, Gen(P) = G1(P) U Supp(P). We test at two clauses/atoms 
ratios. The first test is at 4.258 which is in the phase transition region 
where roughly 50% of the generated 3-SAT clause sets are satisfiable. The 
second test is at clauses/atoms ratio 3.750 where practically all generated 
3-SAT clause sets are satisfiable. 

The test results are shown in Figure [2j In the first test set the key 
problem seems to be finding at least one model candidate. The simpler 
generator (GnTI) appears to perform relatively well except for a few outliers, 
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i.e. instances with significantly higher running time than the average. The 
outliers occur when the generator program G1(P) allows a high number of 
candidate models. At clauses/atoms ratio 3.750 the frequency of outliers for 
GnTI increases and outliers occur already in smaller problem sizes. The 
more involved generating program Gen(P) behaves in a much more robust 
way and the average running time of GnT2 is significantly lower than that 
of GnTI. Next we use the same two test sets for comparing GnT2 and 
dlv. The results are shown in Figure 03 The systems scale very similarly 
in both test sets but DLV seems to be roughly a constant factor faster than 
GnT2. This is probably due to the overhead caused by the more complicated 
generating program in GnT2 and by the two level architecture of GnT2 
where two instances of smodels are cooperating. 

5.3 QUANTIFIED BOOLEAN FORMULAS 

We continue the comparison of GnT2 and dlv using instances of quantified 
Boolean formulas (QBFs) and develop a new way to encode such formulas 
as disjunctive logic programs. In our experiments, we consider a specific 
subclass of QBFs, namely 2, 3-QBFs. Such formulas are of the form 3XVY0 
where X and Y are sets of existentially and universally quantified proposi- 
tional variables, respectively, and is a Boolean formula based on X UY . 
Deciding the validity of such a formula forms a S^-complete decision problem 
|T7] even if is assumed to be a Boolean formula in 3DNF Recall that 
checking the existence of a stable model for a disjunctive logic program is of 
equal computational complexity |12| . which implies the existence of polyno- 
mial time transformations between the decision problems mentioned above. 
In fact, Eiter and Gottlob [12] show how a QBF of the form 3XVY<p with 
cf) in 3DNF can be translated into a disjunctive logic program P such that 
3XWY(p is valid if and only if P has a stable model. This translation is used 
by Leone et al. [2H] to compare dlv and GnT2. 

However, we present an alternative transformation in order to obtain a 
better performance for the two systems under comparison. Our transforma- 
tion is based on the following ideas. The first observation is that we can 
rewrite with 4> in DNF as 3X-3Y-i<j) where ->0 can be understood 

as a Boolean formula in CNF or as a set of clauses S so that each clause 
c G S can represented disjunction of the form 

X 1 V ^X 2 V Y x V ~Y 2 (11) 



35 



Disjunctive 3-SAT: clauses/atoms = 4.258 



le+06 F 




100 150 200 250 300 



Atoms 



Disjunctive 3-SAT: clauses/ atoms = 3.750 




Atoms 



Figure 2: Minimal Models: GnTI vs. GnT2 
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Figure 3: Minimal Models: dlv vs. GnT2 
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where X\ and ~X 2 are the sets of positive and negative literals, respectively, 
which appear in c and involve variables from X while Yy and -Y 2 are similarly 
related to Y. It follows that BXVYcj) is valid if and only if we can find an 
interpretation 3 I : XUY — ► {t,f} such that -xfrnx) is unsatisfiable where 
^(fiiix) denotes the set of clauses Y\ V -Y 2 for which ([TTf belongs to S and 
7 ^ Xi V ~X 2 . The second idea behind our transformation is to choose the 
truth value of the condition Xj V -X 2 for each clause ([TIT) in S rather than 
the truth values of the variables in X U Y. This line of thinking leads to the 
following translation of S. 

Definition 5.1 A clause c of the form fTl\) where Xi, X 2 C X and Y±, Y 2 C 

Y is translated into following sets of rules: 



where c and c are new atoms associated with the clause c, and f and u are 
new atoms. A set of clauses S is translated into 



We use Boolean variables from X U Y as propositional atoms in the trans- 
lation. Intuitively, the rules of Try(c) choose whether a clause c is active, i.e. 
X\ V -1X2 evaluates to false so that the satisfaction of the clause c depends on 
the values assigned to YyUY 2 . The rules in TrE(c) try to explain the preceding 
choice by checking that the values of the variables in X can be assigned ac- 
cordingly. Finally, the rules in Tru(c) implement the test for unsatisfiability 
together with the rule u <— ~w. Basically, the same unsatisfiability check is 
used in the translation proposed by Eiter and Gottlob. However, the trans- 
formation given in Definition 15.11 is more economical as it uses far less new 
atoms and disjunctive rules. In particular, note that variables from X UY 
not appearing in the clauses do not contribute any rules to the translation. 

Next we address the correctness of our transformation and consider a 
2, 3-QBF 3XVV0 where <ft is i n DNF and the disjunctive logic program P 
obtained by translating ->cj) (a set of clauses S) according to Definition 15. 11 

3 The values assigned by / to the variables in Y are not important, but make I a proper 
interpretation over XUY. 
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Lemma 5.2 Let M C Hb(P) 5e a £oia/ propositional interpretation for the 
translation P such that for every clause c G S of the form ( f77|) . (i) c G M 
cgM and (ii) c G M X 1 n M = and X 2 C M. 

Then the programs Try(c), Tte(c), and Tru(c) associated with a clause 
c G S of the form satisfy the following. 

(Rl) The fact c <— belongs to Try(c) M c G M . 

(R2) The fact c <— belongs to Try(c) M -<=>• c G M. 

(R3) For x E Xi, the rule f <— x belongs to Tte(c) A/ -<=>- 
x^M and f G" M. 

(R4) For x G X 2 , the fact x <— belongs to Tte(c) M x G M. 

T/ie ™/e / «- X 2 6elon#s fo Tr E (c) M X 2 M and f £ M. 

(R6) For y G Y"i U 5^, i/ie ra/e y <— u belongs to Tru(c) A/ unconditionally. 

(Rl) The rule Y x U {u} <- Y" 2 6eIon 3 s to T ru (c) M c G M. 

Theorem 5.3 The quantified Boolean formula BXVYcj) is valid if and only 
if the translation P has a stable model. 

PROOF. We may safely assume that all variables inXUF actually appear 
in 0, since redundant variables can be dropped without affecting the validity 
of the formula nor the structure of its translation. 

( => ) Suppose that 3XVF0 is valid. Then there is an interpretation I : 
X U Y -> {t, f } such that I |= VF0. Then define Xj = {x £ X \ I(x) = t}. 
Without loss of generality we may assume that Xj is minimal, i.e. there is 
no interpretation J such that J \= MYcp and Xj C Xj. Then define a total 
propositional interpretation 

M = U Y U {u} U 

{c I c = Xi V -nX 2 V Yi V e ->0 and J ^ X x V -nX 2 } U 

{c I c = Xi V ^X 2 V Y x V ^Y" 2 G ^0 and / |= X x V ^X 2 }. (12) 

It is verified next that M is a stable model of P. The definition of M implies 
that M satisfies the requirements of Lemma l5^2l Then (R1)-(R7) effectively 
describe the structure of P M and it is easy to verify that M is a model of P M 
on the basis of these relationships, as YU {u} C M and / G" M by definition. 
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Next we assume that X C M is a model of P and show that M C X. 
(i) If c G M for some clause c of the form ifTTj) . then c G X, as c <— belongs 
to P M by (Rl). (ii) Similarly, c G M implies c 6 iV by (R2). (iii) We have 
u G X because otherwise X would form a model of by satisfying 

the rules Y x U {u} <- F 2 included in P M by (R7). (iv) Moreover, F C X 
holds, as u G X and X satisfies all the rules y <— u belonging to P M by 
(R6). (v) Let us define an interpretation J : X U Y — > {t,f} such that for 
x' G X, J(x') = t x' G iV, and for y G F, J(y) = I(y). Using (R4), 

we can establish for any (fTT| that I \t= Xi V -X 2 implies J ty= X\ V -X 2 . 
Thus -«f) i(x) ^= ~ [ ^>j{X) where ~^4>i(x) is known to be unsatisfiable. The same 
follows for -i0 j(x) so that J qualifies as an assignment for which J \= VF0 
holds. But then the minimality of I implies J — I, Xj = Xj, and Xj C N. 
To conclude the preceding analysis, M C N and M is a stable model of P. 

( <= ) Suppose that P has a stable model M. Then define an interpreta- 
tion / : X U Y {t, f } by setting = t z G M for any z G X U F. 
Let us then establish that M and / satisfy (fT2ll . (i) The definition of / im- 
plies that Xi = Mfll. (ii) Now u G M, because P contains u <— ~m and 
M is a stable model of P. (iii) For the same reason, / ^ M, because all the 
rules having / as the head have ~/ among the negative body literals, (iv) 
Since u G M and P M contains the rule y <— u for every y G Y, we obtain 
F C M. (v) For any clause c of the form (fTT| . the structure of Try(c) C P 
implies that c G M <^=^> c <— belongs to Try(c) M C P M <^=^> c G" M. 
Using this property, we can establish that c G M <^=^ / ^ Xi V ~X 2 holds 
for the interpretation / defined above, (vi) Thus c G M <^=^ / |= X x V -iX 2 
is implied by the fact that c G M <^=^> c G^ M, as shown above in (v). 

It remains to show that is unsatisfiable. So let us assume the 

contrary, i.e. there is a model Y' C F for -«f>i(x)- Note that M meets the 
requirements of Lemma, l5.2l bv (v) and (vi) above, as the definition of / implies 
iy=X x V ^X 2 X x n M = and X 2 C M. The relationships (R1)-(R7) 

imply that X = (M - (Y\j{u}))UY' is a model of P M , too. Since u £ X, we 
have N C M indicating that M is not stable, a contradiction. Thus -^4>i(x) 
is unsatisfiable which implies I \= VY(f> and the validity of 3XVF0. □ 

Using an implementation of the translation given in Definition 15.11 we 
are able to transform 2, 3-QBFs into disjunctive programs. The remaining 
question is how to generate 2, 3-QBF instances. We use two different schemes 
based on random instances jHllHl. In the first scheme, the sets of variables 
X and F satisfy |X| = |F|. Each random instance is based on v = \X\ + 
\Y\ variables and a Boolean formula which is a disjunction of d = 2 x v 
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conjunctions of 5 random literals out of which at least two literals involve 
a variable from Y, as suggested by Gent and Walsh jTH]. This scheme is 
slightly different from 2QBF GW in [25j based on 3 literal conjunctions just to 
obtain a more challenging benchmark. The constant factor 2 in the equation 
relating d and v has been determined as a phase transition point for the dlv 
system by keeping v — 50 fixed and varying the number of disjunctions in </>. 
In the actual experiment, the number of v variables is varied from 5 to 50 
by increments of 5. We generate 100 instances of 2, 3-QBFs for each value 
of v and translate them into corresponding disjunctive logic programs. The 
running times for dlv and GnT2 are depicted in the upper graph of Figure 
IU The systems scale very similarly, but DLV is on the average from one to 
two decades faster than GnT2. 

In the second experiment with 2, 3-QBFs, we use a different scheme for the 
number of disjunctions d = [-y/v /2J as well as the number of literals which is 3 
in each conjunction. The resulting instances are much easier to solve, because 
d remains relatively low (e.g. d ~ 41 for v = 3500) and many variables do 
not appear in at all. We let v vary from 50 to 3550 by increments of 50 
and generate 100 instances of 2, 3-QBFs for each value of v. The resulting 
running times are shown in the lower graph of Figure 0] The shapes of the 
curves are basically the same, but the performance of GnT2 degrades faster 
than that of DLV. However, the benefits of the translation given in Definition 
15.11 are clear, as GnT2 is able to solve much larger instances than reported 
in where 40 variables turn out to be too much for GnT2. As far as 
we understand, this is due to the sizes of search spaces associated with the 
translated instances of 2, 3-QBFs. For the translation given in Definition [Hill 
the size of the search space examined by GnT2 is of order 2V / ^ whereas it 
is of order 2 V if the translation proposed by Eiter and Gottlob [TJ] is used. 

5.4 PLANNING 

In order to get an idea of the overhead of GnT2 when compared to SMOD- 
ELS, we study three blocks world planning problems encoded as normal pro- 
grams 

• large. c is a 15 blocks problem requiring a 8 step plan using the encoding 
given in [£4] allowing parallel execution of operators, 

• large. d is a 17 blocks problem with a 9 step plan and 
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Quantified Boolean Formulas: Scheme 1 
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Figure 4: Quantified Boolean Formulas: dlv vs. GnT2 
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Table 1: Planning: SMQDELS vs. GnT2 



Problem 


Number of 


Number of 


Time (s) 


Time (s) 




steps 


ground rules 


SMODELS 


GnT2 


large, c 


8 


81681 


4.5 


10.3 




7 


72527 


0.6 


2.1 


large, d 


9 


127999 


10.1 


21.2 




8 


115109 


1.4 


5.2 


large, e 


10 


191621 


18.2 


35.0 




9 


174099 


2.2 


8.7 



• large. e is a 19 blocks problem with a 10 step plan. 

Tabled] contains two entries for each problem: one reporting the time needed 
to find a valid plan with the "optimal" number of steps given as input and one 
reporting the time needed to show optimality, i.e., that no plan (no stable 
model) exists when the number of situations is decreased by one. The times 
reported for each test case are the execution times of SMODELS and GnT2 
given a ground normal program (generated by LPARSE) as input. The results 
show that there is some overhead in the current implementation of GnT2 
even for normal programs and GnT2 handles these examples 2-4 times slower 
than SMODELS. 

6 CONCLUSIONS 

The paper presents an approach to implementing partial and disjunctive 
stable models using an implementation of stable models for disjunction-free 
programs as the core inference engine. The approach is based on unfolding 
partiality and disjunctions from a logic program in two separate steps. In the 
first step partial stable models of disjunctive programs are captured by total 
stable models using a simple linear program transformation. Thus, reasoning 
tasks concerning partial models can be solved using an implementation of 
total models such as the dlv system. This also sheds new light on the 
relationship between partial and total stable models by establishing a close 
correspondence. In the second step a generate and test approach is developed 
for computing total stable models of disjunctive programs using a core engine 
capable of computing stable models of normal programs. We have developed 
an implementation of the approach using SMODELS as the core engine. The 
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extension is fairly simple consisting of a few hundred lines of code. The 
approach turns out to be competitive even against a state-of-the-art system 
for disjunctive programs. The efficiency of the approach comes partly from 
the fact that normal programs can capture essential properties of disjunctive 
stable models that help with decreasing the computational complexity of 
the generate and test phases in the approach. However, a major part of 
the success can be accounted for by the efficiency of the core engine. This 
suggests that more efforts should be spent in developing efficient core engines. 
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